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Abstract 



We study the production of longitudinally polarized A-baryons in single-spin pp collisions 
at RHIC and HERA-iV as a means of determining the spin-dependent A fragmentation 
functions. It is shown that a measurement of the rapidity distribution of the A's would 
provide an excellent way of clearly discriminating between various recently suggested sets 
of polarized A fragmentation functions that are all compatible with present e + e~ data. 
We also address the main theoretical uncertainties, which appear to be well under control. 



The understanding of spin-dependent deep-inelastic scattering processes (DIS) in terms 
of QCD-evolved polarized parton distributions Af(x,Q 2 ) (/ = q,q,g) is still far from 
being satisfactory, despite significant experimental and theoretical progress over the past 
few years. In particular the angular momentum component of the proton's spin and the 
polarized gluon density Ag(x, Q 2 ) remain almost completely unknown for the time being, 
and more experimental results are required. 

Studies of spin-transfer reactions could provide further invaluable and completely new 
insight into the field of "spin physics" and, in addition, might also yield a better under- 
standing of the hadronization process. Such cross sections can be expressed as convolu- 
tions of perturbatively calculable partonic spin-transfer cross sections with certain sets of 
parton distributions and fragmentation functions, whose scale dependence is completely 
predicted by QCD once a suitable non-perturbative input at some reference scale has been 
determined by data. To obtain a non- vanishing twist- 2 spin-transfer asymmetry, the mea- 
surement of the polarization of one outgoing particle is obviously required, in addition to 
having a polarized beam or target. This certainly provides a great experimental challenge. 
A-baryons are particularly suited for such studies due to the self-analyzing properties of 
their dominant weak decay A — ► p7r~ , and recent results on A production reported from 
LEP PJ have demonstrated the experimental feasibility of successfully reconstructing the 
A spin. 

In [[| a first attempt was made to determine the spin- dependent A fragmentation func- 
tions by analyzing these LEP data |l| in leading and next-to-leading order QCD, using 
the results of a preceding study of unpolarized A fragmentation functions. Unfortunately 
it turned out, however, that the available LEP data, all obtained on the Z resonance and 
hence only sensitive to the flavor non-singlet part of the cross section, cannot even suffi- 
ciently constrain the valence fragmentation functions for all flavors. Rather different, but 
all physically conceivable, scenarios adopted for the input valence distributions appear to 
describe the data equally well, and for the "unfavoured" sea quark and gluon fragmenta- 
tion functions one has to fully rely on mere assumptions. Clearly, further measurements 
of other helicity transfer processes are required to test the models proposed in [[J . 
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With the advent of RHIC j3], spin transfer reactions can be studied for the first 
time also in pp scattering at c.m.s. energies of up to ^/s = 500 GeV. In the following 
we will demonstrate that such measurements would provide a particularly clean way of 
discriminating between the various conceivable sets of spin-dependent A fragmentation 
functions presented in 0, and are almost unaffected by theoretical uncertainties. For this 
purpose, only one polarized beam at RHIC would be needed. It should be noted here that 
similar (and almost equally useful) measurements could be performed also in a possibly 
forthcoming experiment at DESY, HERA-iV [[|, utilizing the existing polarized "fixed" 
gas target of HERMES and the impolarized HERA proton beam. 

The process we are interested in is pp — > AX, the arrows denoting a longitudinally 
polarized particle. For the time being, the required partonic helicity transfer cross sec- 
tions, i.e., qq — > qq, . . ., gg —> gg, are calculated only to leading order (LO) accuracy and 
can be found, for example, in |J. Hence we have to restrict our analysis to LO, implying 
the use of LO-evolved A fragmentation functions, contrary to the case of e + e~ — ► AX or 
SIDIS (ep — > AX, ep — > AX) where all relevant coefficient functions are now available 
at next-to- leading order (NLO) || [7[ [|, 0. In various analyses of processes sensitive to 
polarized parton distributions it has turned out to be particularly useful to study distri- 
butions differential in the rapidity of a produced particle ||, to which we therefore limit 
ourselves also in the present analysis. 

The relevant differential polarized cross section can be schematically written as (the 
subscripts "+" , "— " below denote helicities) 

dAa pp ~^ Kx da pp +^ A+x da pp ~^ A + x 



drj drj drj 

dAa ff'^x 
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the sum running over all possible LO subprocesses, and where we have integrated over 
the transverse momentum px of the A, with p™ m denoting some suitable lower cut-off to 
be specified below. The Af p (f p ) are the usual (un)polarized parton distributions of the 
proton, and 

AA A (^ 2 ) = D?$(z,^) - Dfcfaif) (2) 
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describe the fragmentation of a longitudinally polarized parton i into a longitudinally 
polarized A, where Df/^?(z : /i 2 ) (DfS~?(z, jj 2 )) is the probability for finding a A at a mass 
scale fi with positive (negative) helicity in a parton i with positive helicity, carrying a 
fraction z of the parent parton's momentum. 

The directly observable quantity will be not the cross section in ([!]) itself but the 
corresponding spin asymmetry, defined as usual by 

- daPP^ AX /drj { } 

where the unpolarized cross section da pp ^ AX / drj is given by an expression similar to the 
one in ([I]), with all A's removed. 

To study the sensitivity of ([|) to the poorly known A fragmentation functions AD A , 
we use the three LO sets obtained in 0. For the discussion below the idea behind these 
very different models for spin-dependent A fragmentation should be briefly recalled here: 

Scenario 1 is based on expectations from the non-relativistic naive quark model, where 
only strange quarks can contribute to the fragmentation processes that eventually yield 
a polarized A. 



Scenario 2 is inspired by estimates of Burkardt and Jaffe [|1(], [TT[ for a fictitious DIS 
structure function g A , taking into account a similar breaking of the Gourdin-Ellis- Jaffe 



sum rule |T2| for A's as is observed for nucleons. Assuming the same features also for 
the AD A , a sizeable negative contribution from u and d quarks to A fragmentation is 
predicted here. 

Scenario 3 is the most extreme counterpart of scenario 1 since all the polarized fragmen- 
tation functions are assumed to be equal here, which might be realistic if, for instance, 
a sizeable contribution to the production of polarized A's results from decays of heavier 
hyperons who have inherited the polarization of u and d quarks produced originally. 

For the unpolarized parton distributions of the proton, f p , appearing in (|ip, (Q) we use 
the LO set of Ref. [pi throughout our calculations (using other recent LO sets would not 



3 



lead to any sizeable differences here). Unless otherwise stated we use for the corresponding 
polarized densities Af p the LO GRSV "standard" scenario ||14|| . For the unpolarized 
A fragmentation functions D A needed for calculating da pp ^ AX / dr] we use the LO set 
presented in ||, which provides an excellent description of all available, rather precise 
e + e _ data. It should be emphasized, however, that there are still sizeable uncertainties 
for the D A , mainly related to possible SU(3)f breaking effects not discernible from the 
presently available data. We note that in contrast to this, the assumption of SU(2)f 
symmetry (D A = D A ) appears to have a far more solid foundation. Clearly, further 
measurements of the D A are required here. Nevertheless, the uncertainty in the D A 
resulting from SU (3) / breaking does not really affect our conclusions to be drawn below, 
since the contribution from strange quark fragmentation to the unpolarized cross section 
is only about 5%. 

Fig. 1(a) shows our predictions for the spin asymmetry A A as a function of rapidity, 
calculated according to Eqs. (g) and (|TJ) for y/s = 500 GeV and p^ in = 13 GeV. Note that 
we have counted positive rapidity in the forward region of the polarized proton. We have 
used the three different scenarios for the AD A discussed above, employing the hard scale 
fi = Pt- The "error bars" should give an impression of the achievable statistical accuracy 
for such a measurement at RHIC. They have been estimated via 

A4 A ~1 — 1 , (4) 

assuming a polarization P of the proton beam of about 70%, a branching ratio b\ = 
BR(A — > pn) ~ 0.64, a conservative value for the A detection efficiency of e\ = 0.1, and 
an integrated luminosity of C = 800 pb -1 [JJ]. The cross section a pp ^ AX is the unpolarized 
one, integrated over suitable bins of rj. It should be mentioned that results almost identical 
to the ones in Fig. 1(a) can be obtained also for a lower c.m.s. energy of ^/s = 200 GeV 
and a correspondingly lowered p^ tn and luminosity of 8 GeV and 240 pb -1 , respectively. 
Fig. 2(a) shows our results for a conceivable future measurement at HERA-A^ at a much 
lower energy yfs = 40 GeV and for p^ in = 4 GeV and C = 240 pb -1 0. It should be 
stressed that the pt cuts we have introduced do not only guarantee the applicability 
of perturbative QCD (the hard scale fi in (P should be 0(pt)), but also ensure that 
finite-mass corrections to the cross section, which would become increasingly important 



for z < 0.05, remain small [Q]. Furthermore, small values of z also have to be excluded 
in order to make sure that there are no unreasonably large NLO contributions: as was 
noticed for the DIS case in ||, the (unpolarized) NLO kernels for the evolution of the 
fragmentation functions have an extremely singular behaviour at small z, which eventually 
must lead to a complete breakdown of the "perturbative" formalism we use. 

The behaviour of A A in Figs. 1(a) and 2(a) for the different sets of polarized A fragmen- 
tation functions can be easily understood from the fact that the process, in this particular 
kinematical region, is dominated by contributions from u and d quarks, so that the differ- 
ences between the predictions in Figs. 1(a) and 2(a) are driven by the differences in the 
corresponding AD A and AD§. This immediately implies that the asymmetry has to be 
close to zero for scenario 1, negative for scenario 2 and positive and larger for scenario 3. 
The ^-dependence is also readily understood: at negative rj, the parton densities of the 
polarized proton are probed at small values of x<i (i.e., in the "sea region"), where the 
ratio Aq(x 2 ) / q(x 2 ) is also small. On the contrary, at large positive rj, typical values of 
X2 correspond to the valence region where the quarks are polarized much more strongly, 
resulting in an asymmetry that increases with rj. 

The results in Figs. 1(a) and 2(a) clearly demonstrate the usefulness of such kind of 
measurements to determine the polarized A fragmentation functions more precisely. The 
expected statistical errors are much smaller than the differences in A A induced by the 
various models. Thus an analysis of A A would provide an excellent way of ruling out some 
of the presently allowed sets of spin-dependent A fragmentation functions, provided the 
observed differences in A A are not obscured or washed out by the theoretical uncertainties 
inherent in this calculation. We will therefore finally address this important point in 
some detail to demonstrate that the uncertainties appear to be well under control for this 
particular process and do not impose any severe limitations. 

There are three major sources of uncertainties: the dependence of A A on variations 
of the hard scale \x in (|I]), which is of particular importance since we are limited to a 
LO calculation, our present inaccurate knowledge of the precise x-shape and the flavor 
decomposition of the polarized densities Af p , especially of Ag, and our ignorance of 
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AD A which is not constrained at all by the presently available e + e data [fj]. Fig. 1(b) 



gives an example of the scale dependence of A A by changing the scale from /x = p? to 
(jl = Pt/2 for scenario 3. The same is shown in Fig. 2(b) for the HERA-iV situation and 
scenario 2. Even though dAa/drj and da/dt] individually change by as much as a factor 
2 at certain values of rj, the uncertainty almost cancels in the ratio A A . This gives us 
some confidence that the (unknown) NLO corrections might also cancel to some extent in 
the asymmetry, a pattern observed for all available NLO corrections involving polarized 
particles. We also show in the same figures the changes in the predictions resulting from 
varying the polarized parton distributions, using the recent LO set 1 of Ref. |15||, denoted 



by DSS, instead of the GRSV (JT4J one. As can be observed, the asymmetry remains 
practically unchanged, and differences can only be found at the end of phase space (at 
large values of rj) where the cross section becomes small anyway. Also, as an extreme way 
of estimating the impact of the polarized gluon distribution, we have artificially set it to 
zero (Ag(x,fi 2 ) = 0). We find that changes in our predictions only occur in the region 
of negative rj, but are small in the interesting region rj > where the asymmetries are 
larger. 

Finally, in order to examine the role played by AD A in our analysis, we have used two 
different approaches: the standard one for our polarized fragmentation functions, where 
the polarized gluon fragmentation function is assumed to be vanishing at the initial scale 
and is then built up by evolution ("std. AD^), and a set corresponding to assuming 
AD A = D A at the same initial scale of Ref. || ( "max. AD A " ) while keeping the input 
quark fragmentation functions unchanged. As can be observed, the resulting differences 
are also negligible, again due to the fact that u and d fragmentation dominate. 
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Figure Captions 

Fig. 1 (a) The asymmetry A A as defined in (|3|) as a function of rapidity of the A at 
RHIC energies for the various sets of spin-dependent fragmentation functions. The 
error bars have been calculated according to (^) and as discussed in the text, (b) 
same as for scenario 3 in (a), but using the "maximal" AZ}^ (see text), a hard scale 
fi = Pt/2, = 0, or the spin-dependent parton distributions of the proton of set 
1 of fT5| . For comparison the solid line repeats the original result for scenario 3 of 
(a). 

Fig. 2 (a) Same as Fig. 1(a), but for HERA- A kinematics, (b) Same as Fig. 1(b), but 
for HERA- A kinematics and scenario 2. 
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